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This paper discusses a generalization of stimulated Raman adiabatic passage (STIRAP) in which 
the single intermediate state is replaced by N intermediate states. Each of these states is connected 
to the initial state ipi with a coupling proportional to the pump pulse and to the final state ipf with 
a coupling proportional to the Stokes pulse, thus forming a parallel multi-A system. It is shown that 
the dark (trapped) state exists only when the ratio between each pump coupling and the respective 
Stokes coupling is the same for all intermediate states. We derive the conditions for existence of a 
more general adiabatic-transfer state which includes transient contributions from the intermediate 
states but still transfers the population from state ipi to state ipf in the adiabatic limit. We present 
various numerical examples for success and failure of multi-A STIRAP which illustrate the analytic 
predictions. Our results suggest that in the general case of arbitrary couplings, it is most appropriate 
to tune the pump and Stokes lasers either just below or just above all intermediate states. 
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I. INTRODUCTION 

Stimulated Raman adiabatic passage (STIRAP) @-§] 
is an established technique for efficient population trans- 
fer in three-state systems in A or ladder configurations. 
In the original STIRAP, the population is transferred 
adiabatically from an initial state ipi to a final target 
state ipf via an intermediate state ipint by means of two 
partly overlapping laser pulses, a pump pulse flp(t) link- 
ing states ipi and ipi n t, and a Stokes pulse fls(t) linking 
states ipmt and ipf. By applying the Stokes pulse be- 
fore the pump pulse (counterintuitive pulse order) and 
maintaining adiabatic-evolution conditions, one ensures 
population transfer from the initial state into the final 
state, with negligible population in the intermediate state 
at any time. This is so because the transfer is real- 
ized via an adiabatic state ipo(t) (instantaneous eigen- 
state of the Hamiltonian) — the so-called dark state — 
which is a linear superposition of states ipi and ipf only. 
In the ideal limit, unit transfer efficiency is guaranteed 
and the process is robust against moderate changes in 
the laser parameters. Various aspects of STIRAP have 
been studied in detail theoretically and experimentally 
H . Among them are the effects of intermediate-state de- 
tuning Q and loss rate |0L two-photon detuning 
nonadiabatic effects [[l0|-fi2[, multiple intermediate and 
final states |l3)-|l6[] . The success of STIRAP has encour- 
aged its extensions in various directions, such as popula- 
tion transfer in chainwise connected multistate systems 
fl7| prjf and population transfer via continuum 27 - 35 1 . 

In the present paper, we examine the possibilities to 
achieve complete adiabatic population transfer in the 
case when the single intermediate state in STIRAP is re- 
placed by N states, each of which is coupled to the initial 
state ipi with a coupling proportional to the pump field 
and to the final state ipf with a coupling proportional to 



the Stokes field, thus forming a parallel multi-A system. 
In the first place, our work is motivated by the possi- 
bility that appreciable single-photon couplings to more 
than one intermediate states can exist in a realistic phys- 
ical situation, for example, when the pump and Stokes 
lasers are tuned to a highly excited state in an atom, or 
in the case of population transfer in molecules. Such cou- 
plings may be present because, while very sensitive to the 
two-photon resonance @,|, STIRAP is relatively insensi- 
tive to the single-photon detuning from the intermediate 
state j^] , the detuning tolerance range being proportional 
to the squared pump and Stokes Rabi frequencies. It is 
therefore important to know if STIRAP-like transfer can 
take place in such systems. A second example for multi-A 
systems can be found in population transfer in multistate 
chains. It has been shown that when the couplings be- 
tween the intermediate states arc constant the multistate 
chain is mathematically equivalent to a multi-A system, 
in which the initial state ipi is coupled simultaneously 
to N — 2 dressed states which are in turn coupled to 
the final state ipf |2(|. It has been demonstrated that 
in certain domains of interaction parameters (Rabi fre- 
quencies and detunings), adiabatic population transfer 
in these multistate chains (respectively, in the equivalent 
multi-A systems) can take place, while in other domains 
it cannot. A third motivation for the present work is 
population transfer via continuum |2^-|35||. In their pio- 
neering work on this process p7| , Carroll and Hioe have 
replaced the single descrete intermediate state in STI- 
RAP by a quasicontinuum consisting of infinite number 
of equidistant discrete states with energies going from 
—00 to +00. Moreover, each of these states was coupled 
with the same coupling Qp(t) to the initial state ipi and 
with the same coupling fis(t) to the final state ipf. Un- 
der these conditions Carroll and Hioe have shown that 
complete population transfer is achieved in the adiabatic 
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limit with counterintuitively ordered pulses. It has been 
shown later |nj that the Carroll-Hioe's quasicontinuum 
is too simplified and symmetric, and that a real contin- 
uum has properties, such as nonzero Fano parameter, 
which prevent complete population transfer. It is inter- 
esting to find out which of the simplifying assumptions in 
the Carroll-Hioc model (equidistant states, going to infin- 
ity both upwards and downwards, equal pump couplings 
and equal Stokes couplings) makes the dark state in the 
original STIRAP remain as a zero-eigenvalue eigenstate 
of the multi-A Hamiltonian. In this paper, we consider 
the general asymmetric case of unequal couplings and 
unevenly distributed, finite number intermediate states. 
Besides the Carroll- Hioe model [27|] , our paper general- 
izes the results of Coulston and Bcrgmann 131 who were 
the first to consider the effects of multiple intermediate 
states in the simplest case of N = 2 states and equal 
couplings Qp(t) to state ipi and equal couplings ils(t) 
to state ipf. We shall find the conditions for the exis- 
tence of the dark state in multi-A systems as well as the 
conditions for the existence of a more general adiabatic- 
transfer state, which still links states ipi and ipf adiabat- 
ically but is allowed to contain transient contributions 
from the intermediate states. 

Our paper is organized as follows. In Sec. || we present 
the basic equations and definitions and review the stan- 
dard three-state STIRAP. In Sec. [nj we discuss the case 
when all intermediate states are off single-photon reso- 
nance. In Sec. IV we consider the case when one of the 



intermediate states is on single-photon resonance and in 
Sec. [v] the case of degenerate resonant states. In Sec. Wj\ 
we use the adiabatic-elimination approximation to gain 
further insight of the process. Finally, in Sec. VII 
summarize the conclusions. 
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II. BASIC EQUATIONS AND DEFINITIONS 

A. Basic STIRAP 

The probability amplitudes of the three states in STI- 
RAP satisfy the Schrdinger equation (h = 1), 



i- c (t) = n(t)c(t), 



(i) 



where c(t) — [ci(t),Ci nt (t),Cf(t)] T . In the rotating-wave 
approximation pq] , the Hamiltonian is given by 



H(t) = 








A 



o 

n s (t) 
o 



(2) 



The time- varying Rabi frequencies £lp(t) and £ls(t) are 
given by products of the corresponding transition dipole 
moments and electric- field amplitudes. States ipi and ipf 
are assumed to be on two-photon resonance, while the in- 
termediate state ipint may be off single-photon resonance 



by a detuning A. The system is initially in state ipi, 
Ci{— oo) = 1, Ci nt {— oo) = cj(— oo) = 0, and the quanti- 
ties of interest are the populations, and partucularly the 
population of state ipf at t — > +00, Py(oo) = |c/(oo)| 2 . 

Throughout this paper, we assume that the two pulses 
are ordered counterintuitively, i.e., the Stokes pulse pre- 
cedes the pump pulse, 



n P (t) 
n s (t) 



but we do not impose any restrictions on the particu- 
lar time dependences in our analysis. In the numerical 
examples, we assume Gaussian shapes, 

n P (t) = n Q e-^r/T^ n s(t ) = n e-^/ T2 , (4) 

where T is the pulse width, 2t is the time delay between 
the pulses, and we take r = 0.5T everywhere. Further- 
more, we choose the peak Rabi frequency f2o to define 
the frequency and time scales. 

The essence of STIRAP is explained in terms of the 
so-called dark (or trapped) state Lpp>(t), which is a zero- 
eigenvalue eigenstate of H(t), 



where Q(t) is the mean-square Rabi frequency, 



n(t) = yjn%(t) 



n 2 P {t). 



(5) 



(6) 



For counterintuitively ordered pulses, Eq. (J3J) , we have 
<Pd(— oo) = ipi and (pp>(+oo) = —ipf and hence, the dark 
state connects adiabatically states ipi and ipf. By main- 
taining adiabatic evolution (a condition which amounts 
to requiring that the pulse width T is large or that the 
pulse areas are much larger than 7r), one can force the 
system to remain in the dark state and achieve complete 
population transfer from ipi to ipf . Moreover, since y>D{t) 
does not involve the intermediate state ipint , the latter is 
not populated in the adiabatic limit, even transiently, 
and hence, its properties, including decay, do not affect 
the transfer efficiency. 



B. Multi-A STIRAP 

1. The system 

In the multi-A generalization of STIRAP the single 
intermediate state is replaced by N intermediate states 
ipi, ip2, ■ ■ ■ , ipN, each of which is coupled to both states ipi 
and ipf, as shown in Fig. [l| The system is again assumed 
to be initially in state ipi, 



Ci(-00) = 1, 

c/(— oo) = ci(— oo) = 



= ctv(-oo) = 0, 



(7a) 
(7b) 
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FIG. 1. The multi-A system studied in this paper. The 
system is initially in state ipi and the objective is to transfer it 
to state ipf via one or more intermediate states ipi, ip2, ■ ■ ■ , ipN 
by means of two delayed laser pulses, the pump pulse Qp(t) 
and the Stokes pulse fis(i), ordered counterintuitively (the 
Stokes pulse coming first). Each intermediate state is cou- 
pled to state ipi with a coupling Slp^(t) proportional to Qp(t) 
and to state ipf with a coupling Qs,k(t) proportional to Qs(t)- 
Each intermediate state tpt is detuned from single-photon res- 
onance with a detuning At, while states ipi and ipf are on 
two-photon resonance. 

and the objective is to transfer the population to state 
ipf. The column vector of the probability amplitudes 
in the Schrdinger equation (Q) is given by c(i) = 
[ci(t), ci(t), . . . , cat(£), Cf(t)) T and the Hamiltonian reads 



H 
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Qp i2 A 2 
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(8) 



States ipi and ip / are again assumed to be on two-photon 
resonance while each intermediate state ipk may be off 
single-photon resonance by a detuning . The couplings 
Qp,k(t) of the intermediate states ip k to the initial state 
ipi are proportional to the pump field, while the couplings 
Qs,k(t) °f the intermediate states to the final state tpf are 
proportional to the Stokes field, 



Jip,fc(t) = a k sip(t), n Stk (t) = I3 k n s {t). 



(9) 



The dimensionless numbers a k and (3 k characterize the 
relative strengths of the couplings while flp(t) and £ls(t) 
are suitably chosen "units" of pump and Stokes Rabi fre- 
quencies. We fix these "units" by choosing a\ = (3i = 1, 



which means that £lp(t) and $ls(t) axe the pump and 
Stokes Rabi frequencies for state ipi, Qp(t) = ttp } i(t) 
and ils(i) = £ls.i(t), and the constants a k and (3k are 
the relative strengths of the pump and Stokes couplings 
of state ipk with respect to those for state ip\. In real 
physical systems these constants contain Clebsch-Gordan 
coefficients and Franck-Condon factors. Moreover, with- 
out loss of generality we assume that Qp(t), f2s(t), a k , 
and (3k are all positive. 



2. Adiabatic-transfer state 

As has been emphasized in |25|,|26j , a necessary condi- 
tion for complete adiabatic population transfer in mul- 
tistate systems is the existence of an adiabatic-transfer 
(AT) state </St(£) which is defined as a nondegenerate 
eigenstateF] of H(t) having the property 



ip T {t) 



ipi, t — » -oo 
ipf, t — > +oo 



(10) 



up to insignificant phase factors. The dark state (@), 
which is a coherent superposition of states ipi and ipf 
only, is the simplest example of such a state. Under 
certain (quite restrictive) conditions it is an eigenstate 
of the multi-A system too, but in the general case the 
Hamiltonian (|J) does not have such an eigenstate. Un- 
der some more relaxed conditions, however, H(t) has an 
eigenstate with the more general properties ([L0|) (which 
allow for nonzero transient contributions from the inter- 
mediate states), and we derive these conditions below. 



III. THE OFF-RESONANCE CASE 

If all single-photon dctunings are nonzero, A& (k 
1,2,..., N), we have 



det h = n 2 P n 2 s v (s a 2S P 2 - sip) 



(11a) 



N N 



H V kl {a k (3i - ai (3 k ) 2 , (lib) 



fc=i i=k+i 



where 



N , 

at 



N 



5 « 2 =E5;< v=E^. 5 ^=£^> ( 12 ) 

fc=l K k=l K fe=l K 



*If the eigenstate <fiT{t) is degenerate, there will be reso- 
nant population transfer to the other such eigenstate(s) with 
a transition probability sin 2 ^A, where A is the pulse area of 
the nonadiabatic coupling between these eigenstates, with in- 
evitable population loss. A is independent on the pulse width 
T and hence, it does not vanish in the adiabatic limit T — > oo, 
unless the nonadiabatic coupling is identically equal to zero. 
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N 



N 



N 



V=]]_A k , V„=l[A k , V mn = \[ A fc . (13) 



fe=i 



k=i 

k^n 



k = l 



Hence, detH ^ in the general case. This means that, 
unlike in STIRAP (N = 1), H(t) does not necessar- 
ily have a zero eigenvalue. We shall consider first in 
Sec. Ill A the case when a zero eigenvalue exists (with 
the anticipation, in analogy with STIRAP, that the cor- 
responding eigenstate is the desired AT state) and then 
in Sec. IIIB the case when it does not exist. 



A. A zero eigenvalue 

1. Condition for a zero eigenvalue 
The condition for a zero eigenvalue is given by 

N N I a a \2 
b a 2b 2 - h afj = ^ 2^ = °- ( 14 ) 



k=l l=k+l 



A fe A, 



Obviously, this condition depends only on the relative 
coupling strengths and the detunings, but neither on time 
nor on laser intensities. Hence, it remains unchanged as 
the adiabatic limit is approached. 

The eigenstate corresponding to the zero eigenvalue 
most generally reads 



N 



Pa 



(t) = ai{t)ipi + a f (t)tp f + J2 a k(t)^k- (15) 



fc=i 



The amplitudes of the intermediate states are given by 
a>k{t) = r adt) a/(t), (16) 



with k — 1, 2, . . . , N. Obviously, they may be nonzero 
at finite times but vanish at ±oo, afc(±oo) = 0, because 
both £lp t k(t) and fls,k(t) vanish at infinity. The ampli- 
tudes of the initial and final states satisfy both equations 



S a 2n P (t) ai (t) + S al3 n s (t)a f (t) = 0, (17a) 
S a0 n P {t) ai (t) + Sfp(ls(t)af(t) = 0, (17b) 

which are linearly dependent because of Eq. (jl4|) . 

We are going to consider two cases: when each term in 
the sum ( |l4[ ) is zero and when the individual terms may 
be nonzero but the total sum vanishes. 



2. Proportional couplings 



When 



21 

Ih 



a 2 



= 1, 



(18) 



each term in Eq. ( |l4| ) vanishes and the zero eigenvalue 
exists regardless of the detunings. Condition (|18|), which 
is essentially a condition on the transition dipole mo- 
ments, means that for each intermediate state ipk, the 
ratio £lp,k(t) /£ls,k{t) between the couplings to states ipi 
and ipf is the same and does not depend on k. It fol- 
lows from Eq. (jlq) that S a 2 = Sp2 = S a p = S and 
we find from Eq!~@ that ^p^WO) + n s (iW(t) = 0, 
provided S ^ 0. Then it follows from Eq. (n6h that all 
intermediate-state amplitudes are zero, a,k{t) = 0. It is 
easily seen that the zero-eigenvalue eigenstate coincides 
with the dark state (^) and in the adiabatic limit it trans- 
fers the population from state ipi to state ipf, bypassing 
the intermediate states. Hence, when condition (|lj) is 
fulfilled the multi-A system behaves very much like the 
single A system in STIRAP. This confirms and general- 
izes the conclusions in 13 2?| that complete population 
transfer is possible when all ak and (3k are equal. 



3. Arbitrary couplings 

Suppose now that condition (|l^) is not fulfilled while 
condition (|l4|) still holds. This can be achieved by chang- 
ing the two laser frequencies simultaneously while main- 
taining the two-photon resonance, which corresponds to 
adding a common detuning A to all single-photon detun- 
ings, A k -> A fe +A (k = 1, 2, . . . , N). For N intermediate 
states, there are N — 2 values of A for which condition 
( |l4| ) is satisfied. In the zero-eigenvalue eigenstate (|l5| ) 
the amplitudes of the intermediate states rtlq) are gener- 
ally nonzero. However, if 5 Q 2, Sp2, and S a f3 are nonzero 
we have |a,(— oo)| = 1 and |et/(+oo)| = 1, and hence, this 
eigenstate is an AT state, as defined by Eq. (|l0|). 



4- The case of vanishing S a 2, Sa2, and S, 



a/3 



Let us now suppose that some of the sums S a 2 , Sp2 
and S a f3 are equal to zero. If S a 2 ^ and 5^2 = 
(whi ch a lso requires that S a p = 0) then it follows from 
Eq. Ql7ap that a,(i) = 0. Hence, we have |a/(±oo)| = 1 
and <£>o(±oo) = ipf (up to an irrelevant phase factor) and 
there is no AT state. A similar conclusion holds when 
S a 2 — and Sp2 ^ 0: then a/(i) = 0, |aj(±oo)| = 1, and 
V?o(±oo) = ipi. 

The case when S a 2 = Sp2 = S a p — is different be- 
cause then, as can easily be shown, the Hamiltonian has 
two zero eigenvalues. Consequently, there will be reso- 
nant transitions between the corresponding degenerate 
eigenstates, even in the adiabatic limit, which again pro- 
hibit complete STIRAP-like population transfer. 

We shall illustrate this conclusion by calculating the 
final-state population for proportional couplings dis)). 
Due to the degeneracy, we have some freedom in choosing 
the corresponding pair of orthonormal eigenstates. As 

one of them, it is convenient to take state (||), tp^\t) = 
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FIG. 2. The final-state population Pf as a function of 
the pulse width T in the case of TV = 3 intermediate states 
for four combinations of coupling strengths a*, and f3k and 
dimensionless detunings 5k ~ Ak/Qo given in the table. 



ipn(t), because y>£>(— 00) = ipi and (pn(+co) = —ipf- The 
other zero-eigenvalue adiabatic state can be determined 
by Gram-Schmidt ortogonalization and is 



n P (t) . , n s (t) 



n(t) 



■V. 



JV 



n(t) 



fc=l 



(19) 



l + fi 2 (*)Ef=i4/A^ 



-1/2 



, where f2(i) is 



with ^(i) 

given by Eq. (|6|). In the adiabatic limit, the population of 
state ipf, which is equal to the probability of remaining 

in the adiabatic state ¥?q , is given by 



P 



f 



(20) 



where ti(t) = ip 1 ^ (t) • (f^ (t) and tantf(i) = n P (t)/n s {t). 

In conclusion, when H(t) has a zero eigenvalue 
STIRAP-like transfer is possible only if 



S a 2 ^ 0, 



V ^ 0, s a() ^ 0. 



(21) 



5. Examples 

In Fig. ||, the final-state population Pf in the case of 
N = 3 intermediate states is plotted against the pulse 




-3-2-1 12 3 

Time (units of I) 

FIG. 3. The time evolutions of the populations in the case 
of TV = 3 intermediate states for a pulse width T — 30^ \ 
The interaction parameters for the upper and lower plots are 
given respectively in the first and second rows of the table in 
Fig- 1 



width T for four combinations of coupling strengths and 
detunings. The solid curve is for a case when condition 
( |l8| ) is satisfied and the zero-eigenvalue eigenstate is the 
dark state (§). The dotted curve is for a case when con- 
dition (|l8| ) is not satisfied but conditions ( [Til) and (|2l| ) 
are and the zero-eigenvalue eigenstate is an AT state, 
Eq. (|l0|). In both cases, the final-state population Pf 
approaches unity as the pulse width increases and the 
excitation becomes increasingly adiabatic. The dashed 
curve is for a case when S a 2 — Sp2 = S a p = 0; then, as 
the adiabatic limit is approached, Pf tends to the con- 
stant value Pf w 0.442, predicted by Eq. (|o|). Finally, 
the dashed-dotted curve is for a case when Sp2 = <S Q(3 = 
and S a i ^ 0; then, in agreement with our analysis, Pf 
tends to zero in the adiabatic limit. 

In Fig. [j], we show the time evolutions of the pop- 
ulations in the case of TV = 3 intermediate states un- 
der almost adiabatic conditions. The upper plot is for 
the solid-curve case in Fig. when the zero-eigenvalue 
eigenstate is the dark state (|5|); consequently, the inter- 
mediate states remain virtually unpopulated. The lower 
plot is for the dotted-curve case in Fig. ^ when the zero- 
eigenvalue eigenstate is an AT state with nonzero com- 
ponents from the intermediate states; consequently, these 
states acquire some transient populations. 
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B. No zero eigenvalue 

When condition (|l4j) is not satisfied, dct H^O and the 
Hamiltonian (||) does not have a zero eigenvalue. This 
fact alone does not mean much because any chosen eigen- 
value can be made zero by shifting the zero energy level 
with an appropriate time-dependent phase transforma- 
tion. More importantly, an AT state ipx, as defined by 
Eq. ([To]) , may or may not exist and the conditions for its 
existence are derived below. The derivation is similar in 
spirit to that for multistate chains in |^5| , p6| . 

By setting Hp = Hs = in Eq. (||), we find that there 
are two eigenvalues which vanish as t — > ±00 (although 
they are nonzero at finite times), while the others tend 
to the (nonzero) detunings At ±00, each of the two 
eigenstates corresponding to the vanishing eigenvalues is 
equal to either state ipi or state tpf or a superposition of 
them. Obviously, if an AT state exists, its eigenvalue At 
should be one of these eigenvalues. Hence, we have to find 
the asymptotic behaviors of these eigenvalues and the 
corresponding eigenstates, i.e., we need to determine how 
the degeneracy of these eigenvalues is lifted by the laser 
fields. We note here that only one eigenvalue vanishes 
when Hp — * and Hs 7^ 0, which happens at certain 
early times, or when Hs — * and Hp ^ 0, which happens 
at certain late times. 



From here we find Aj'(0), replace it in the Taylor expan- 
sion of A ,7 (f^p), and keeping the lowest-order nonzero 
term only, we obtain 



a: 



s 



a/3 



73= 



Hp. 



(25) 



In order to find the dependence of A ; ~ on H s , we set 
Hp = in Eq. ([22]), divide by A (which amounts to re- 
moving the root Aj), differentiate with respect to f2|, 
set 0| = and \f (Sl 2 s = 0) = 0, and find 



2. Late-time eigenvalues 



(26) 



In a similar way, we find that at late times, when 
Hs/Hp — * 0, the two vanishing eigenvalues behave as 



At 



A^ ~ ^Q,2^1p. 



(27) 
(28) 



1. Early-time eigenvalues 

Let us first consider the case of early times (t — > —00). 
It follows from the above remarks that as soon as the 
pulse Hs arrives, one of the vanishing eigenvalues, A ; ~~ 
(the "large" one), departs from zero, while the other, A J 
(the "small" one) , remains zero until the pulse H p arrives 
later. Since A~ vanishes when Hp — ► and A ; ~~ vanishes 
when Clg — > 0, AJT should be proportional to some power 
of Hp and A ; ~ to some power of Hs ■ Since at those times 
Hp/Hs — * 0, the relation |A~| <C \X~[\ holds; hence, the 
names "small" and "large". 

To determine X~ , we consider the eigenvalue equation, 

det(H - Al) = h + h 1 X+ ...+ h N+2 X N+2 = 0, (22) 

as an implicit definition of the functional dependence of 
Xj on Hp. Note that ho = detH. Since all hk depend 
on flp only via H Pl A J has a Taylor expansion in terms 
of Hp. We differentiate Eq. (p2|) with respect to Hp, set 
Hp = and X~(Hp = 0) = 0, and obtain 

h' o (0) + h 1 (0)X-'(0)=0, (23) 

where a prime denotes d/dHp and 

ti (0) = nlV(S a2 S^-Slp) ) (24a) 
hi(0) = H%VS 2. (24b) 



3. Connectivity and AT condition 

It is easy to verify that the eigenstates corresponding 
to A ,7 and A^ coincide with state tpi, while those corre- 
sponding to A ; ~ and A+ coincide with state tpf. Hence, 
the AT state <px, if it exists, must have an eigenvalue 
that coincides with Aj at early times and with A+ at 
late times. It should be emphasized that A~ and A+ do 
not necessarily correspond to the same eigenvalue and it 
may happen that A~ is linked to A ; + rather than A+; then 
an AT state does not exist. In any case, the upper (the 
lower) of the two eigenvalues at —00 is connected to the 
upper (the lower) of the two eigenvalues at +00. Since 
AH 3> I AT/ 1 and (A^ I » \Xf\, the linkage is determined 
by the signs of the "large" eigenvalues A ; _ and A^~ . If they 
have the same signs, they will be both above (or below) 
A~ and A+ and hence, the desired linkages A ; ~ A ; + 
and Xj «-> A+ will take place. If Xf and Aj* - have op- 
posite signs they cannot be connected because such an 
eigenvalue will cross the one linking A~ and A+, which 
is impossible. Thus, from this analysis and Eqs. fl26| ) 
and (|2^) we conclude that the necessary and sufficient 
condition for existence of an adiabatic-transfer state is 

S a 2S 2 > 0. (29) 

It is easy to sec that condition (Q for existence of a 
zero eigenvalue, along with condition (|2l]) for existence 
of AT state in this case, agree with condition (^9j). In- 
deed, we have S a 2Sp2 = S 2 ^ > 0. Moreover, the zero 
eigenvalue is reproduced correctly by Eqs. ( pB| ) and fl27|). 
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4- The case of vanishing S a 2 and Sg2 



The derivation of the AT condition (|29|) suggests that 
both sums S a 2 and Spz should be nonzero. Let us ex- 
amine the case when one of them is zero, e.g., Sp2 = 0. 
By going through the derivation that lea ds t o Eq. ( p5| ) 
we find that now, as evident from Eq. ( ^4b| ) , we have 
/ii(0) = 0. It follows from Eq. ( p^ ) that now two, rather 
than one, eigenvalues vanish when fip = and Qs ^ at 
early times (because then ho = dctH = too). Hence, 
in contrast to the case of Sp2 ^ 0, the arrival of the 
Stokes pulse Sls(t) does not make one of these eigenval- 
ues depart from zero but rather their degeneracy is lifted 
only with the arrival of the pump pulse ttp(t) later. The 
implication is that the initial state ipi cannot be identi- 
fied with a single adiabatic state at t — > — oo but it is 
rather equal to a superposition of two adiabatic states]]; 
hence, there is no AT state. A similar conclusion holds 
in the case when S a 2 = 0: then the final state ipf cannot 
be identified with a single adiabatic state at t — » +oo. 
Therefore, for S a 2 = or Sa2 = an AT state does not 
exist, as follows formally from Eq. (|2^). 

The situation with the mixed sum S a f3 is different. In 
the above derivation the condition S a p ^ was not re- 
quired anywhere which means that an AT state exists 
even for S a p = 0. We will return to this problem in 
Sec. M 



5. Examples 

In Fig. [| we have plotted the time evolutions of the 
eigenvalues (upper row of figures) in the case of N = 2 in- 
termediate states for two combinations of detunings and 
the same set of coupling strengths. The solid curves are 
calculated numerically and the dashed curves show our 
asymptotic approximations (|25|)-(|2lj|). The bottom row 
of figures show the components of this eigenstate which 
is equal to bare state ipi initially; it corresponds to the 
eigenvalue whose asymptotics at early times is described 
by A~. Hence, the squared components of this eigenstate 
give the populations of the four bare states in the adia- 
batic limit. As we can see in the left column of figures, 
there is an eigenvalue whose asymptotics is given by Aj 
at early times and by A+ at late times; this is so because 
condition ( ^9|) is satisfied in this case. The corresponding 
eigenstate is an AT state, as evident from the bottom left 
figure, because it is equal to state ipi initially and to state 



^The initial state ipi ls associated with a single adiabatic 
state at — oo when the arrival of the pump pulse lifts the 
degeneracy of one and only one eigenvalue. Similarly, the 
final state ipf is associated with a single adiabatic state at 
+oo when the vanishing Stokes pulse restores the degeneracy 
of one and only one eigenvalue. 
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FIG. 4. The upper row of figures show the time evolutions 
of the eigenvalues in the case of N = 2 intermediate states for 
two sets of detunings: Ai = 0.5f2o,A2 = 1.5f2o (left figure) 
and Ai = — 0.5f2o,A2 = 0.5f2o (right figure). The coupling 
strengths are the same in both cases, ot\ = 0i = 1, Q2 = 2, 
and — 0.5. The solid curves are calculated numerically 
and the dashed curves show the asymptotic approximations 
(pB])-(pg|). The bottom row of figures show the components of 
this eigenstate which is equal to the bare state ipi initially. It 
corresponds to the eigenvalue which is shown by a thick curve 
in the corresponding upper figure and whose asymptotics at 
early times is described by AJ. The labels denote the bare 
states to which the components belong and the final-state 
component is shown by a thick curve. 



— ipf finally. For the case shown in the right column of 
figures, there is no AT eigenvalue because the asymptotic 
behaviors Xj and A+ are related to two different eigen- 
values; this is so because condition (|2^ ) is not satisfied 
in this case. Consequently, there is no AT eigenstate, as 
evident from the bottom right figure, because the shown 
eigenstate is equal to state ipi both initially and finally. 

In Fig. H we have plotted the final-state population Pf 
as a function of the pulse width T in the case of N = 2 
intermediate states for the same two sets of interaction 
parameters as in Fig. |. The solid and dashed curves in 
Fig. H correspond to the left and right columns in Fig. |J, 
respectively. As follows from Eq. (p9|), an AT state exists 
for the solid curve and does not exist for the dashed curve. 
Indeed, as seen in the figure, as T increases, the final- 
state population Pf approaches unity for the solid curve 
and zero for the dashed curve. 

In Fig. [6] the final-state population Pf is plotted as a 
function of the single-photon detuning A of the pump and 
Stokes fields from the lowest intermediate state. We have 
taken two intermediate states ipi and ip2 with Ai = A, 
A2 = fig + A. The coupling strengths are taken the 
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FIG. 5. The final-state population P/ as a function of the 
pulse width T in the case of N = 2 intermediate states for the 
same two sets of interaction parameters as in Fig. ^. The solid 
and dashed curves correspond to the left and right columns 
in Fig. ^, respectively. The detunings for the two curves are 
Ai = 0.5flo, A2 = 1.5fio (solid curve), and Ai = — 0.5f2o, 
A2 = 0.5f2o (dashed curve). The coupling strengths are the 
same for both curves, Qi = f3\ = 1, 02 = 2, and fa = 0.5. 
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FIG. 6. The final-state population Pf as a function of the 
single-photon detuning A of the pump and Stokes fields from 
the lowest intermediate state. We have taken two interme- 
diate states ip\ and ^2 with Ai = A, A2 = flo + A. The 
positions of the single-photon resonances with these states 
are shown by the vertical solid lines. The region where an AT 
state does not exist, — 0.8f2o<A< — 0.2f2o, is shown by verti- 
cal dashed lines. The relative coupling strengths are taken as 
in Figs. § and |, ai = fa = 1, a 2 = 2, and fa = 0.5. The thin 
curve is for QoT = 20 and the thick curve is for £IqT = 80. 
The inset shows Pf for QqT — 20 in a wider range. 
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FIG. 7. The final-state population Pf as a function of 
the single-photon detuning A of the pump and Stokes fields 
from the lowest intermediate state. We have taken N = 5 
equidistant intermediate states, i.e., Afe = (k— l)flo + A with 
k = 1, 2, 3, 4, 5. The positions of the single-photon resonances 
with these states are shown by the vertical solid lines. The 
relative coupling strength at and /3k are taken randomly, as 
provided by the random generator of our computer, and given 
in the table (up to four digits after the decimal point). The 
thin curve is for QqT = 20 and the thick curve is for £loT = 80. 
The vertical dashed lines separate the regions where an AT 
state exists or does not exist, as calculated from Eq. (p9[). 

same as in Figs. [| and [| In the region where an AT 
state does not exist, — 0.8J7o < A < — 0.2Hq [calculated 
from Eq. (p9|)], the transfer efficiency is low, while outside 
it the transfer efficiency is almost unity, as a result of the 
existence of an AT state. The left and right columns of 
plots in Fig. |^ (respectively, the solid and dashed curves 
in Fig. ||) correspond to detunings A = 0.5Slo and A = 
— 0.5flo: respectively. The inset shows how the transfer 
efficiency eventually decreases at large detunings which, 
as in STIRAP Q, is due to deteriorating adiabaticity. 

In Fig. 0, we have plotted the final-state population 
Pf as a function of the single-photon detuning A of the 
pump and Stokes fields from the lowest intermediate state 
for the case of N = 5 equidistant intermediate states and 
randomly taken coupling strengths and j3k- As Fig. |?j 
shows, there are three distinct domains of single-photon 
detunings: A < -4O , -4^o < A < 0, and A > 0. For 
A < — 4f2o and A > 0, we have Pf pa 1, whereas for 
— 4f2o < A < 0, there are alternative regions of high and 
low transfer efficiency. This behavior is easily explained 
by the AT condition (|29|). As A changes, we pass through 
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the zero points of the sums S a 2(A) and 5ga(A), thus go- 
ing from an interval where these sums have the same sign 
(where Pt « 1) to an interval where they have opposite 
signs (where Pf ss 0) and vice versa. Obviously, for suf- 
ficiently large and negative A, both <S Q 2(A) and Spz (A) 
are always negative and condition ( p9|) is satisfied, which 
ensures the existence of AT state and STIRAP-like unit 
transfer efficiency. Similarly, for sufficiently large and 
positive A, both <S a 2(A) and Sp2 (A) are always positive 
and we have unit transfer efficiency there too. The con- 
clusion is that whenever the pump and Stokes lasers are 
tuned below or above all intermediate states, STIRAP- 
like transfer is always guarantied in the adiabatic limit. 
When the lasers are tuned within the manifold of in- 
termediate states such a transfer may or may not take 
place, depending on whether the AT condition ( p9| ) is 
satisfied or not. Furthermore, tuning the lasers below 
or above all states appears to be more reasonable even 
than tuning on resonance with an intermediate state be- 
cause adiabaticity is achieved more easily, as the curves 
for Q,qT = 20 in Figs. ||and0 show. Moreover, as follows 
from Eq. (plj|), the transient intermediate-state popula- 
tions P k (t) decrease with the detuning as A^ 2 . 



IV. A RESONANT INTERMEDIATE STATE 



A. A zero eigenvalue 



case of proportional couplings (|l8|), it can readily be 
shown that the zero-eigenvalue adiabatic state is given 
again by the dark state (|J). Moreover, now condition 
( pi| ) is not required and the dark state ipo(t) is a zero- 
eigenvalue eigenstatc of H(t) even when S^' , <Sg™ > ano - 

S^j vanish^. In the case of arbitrary couplings, the sum 
(BOat) can vanish only by accident because we are not al- 
lowed to "scan" the pump and Stokes laser frequencies 
across the intermediate states as we would violate the 
assumed single-photon resonance condition A„ = 0. If 
this happens it can easily be shown that again, the zero- 
eigenvalue eigenstate is an AT st ate, th e only difference 
from the off-resonance case (Sec. Ill A 3| ) being that now 
the AT state does not contain a component from the 
resonant bare state, a n (t) = 0. In both cases, we have 
complete population transfer in the adiabatic limit. 



B. Nonzero eigenvalue 

For dctH 7^ 0, we fo llow the same approach as for 
nonzero detunings (Sec. 1IIB). By setting Up — £1$ — 
in Eq. (||) we find that for A„ = 0, there are three, 
rather than two, vanishing eigenvalues. Hence, we have 
to establish how the new, third, zero eigenvalue affects 
the AT state. 



If a certain detuning is equal to zero, A„ = 0, we have 



detH = n 2 P n 2 s Vnk(a k Pn ~ a n p k ) 2 



k = l 
k^n 



(30a) 



(30b) 



where the iS^ n ^-sums are defined as the 5-sums but with- 
out the n-th terms, 



N 2 

a 



N 



k=l 

k^-n 



k=l 



ft 



N 



a k fik 



k = l 



(31) 



The intermediate-state amplitudes, except a n (t), are 
given by Eq. ( |l6| ) . The equation for a n (t) is replaced by 

«P,„(iH(t) + tt s ,n(t)a f (t) = 0, (32) 

while Eqs. ( |l7| ) are replaced by 

S^np^it) + S { ^n s (t)a f (t) + a n a n (t) = 0, (33a) 

S^n P (t)ai(t) + Spn s (t)a f (t) + (3 n a n {t) = 0. (33b) 

We can again consider the two possibilities: when 
condition (n~|) is fulfilled and when it is not. In the 



1. Early-time eigenvalues 

It is readily seen from Eq. (||) that only one eigenvalue 
vanishes for D,p = and ilg ^ 0. This means that as 
soon as the Stokes pulse f2s(i) arrives, the degeneracy 
of two of the eigenvalues, Aj^ and Xf 2 , is lifted and they 
depart from zero, while the third eigenvalue Aj stays zero 
until the pulse flp(t) arrives later. Hence, there are two 
"large" eigenvalues and one "small" eigenvalue. 

The "small" eigenvalue Aj(fip) can be determined in 
the same manner as for nonzero detunings. We have 

^o(Q) - U? s V n \a 2 n S$ - 2a r ASg + 2 n S%>} , (34a) 



*/3 

hi(Q) = n 2 s plv n . 

Using Eq. (|23|), we find Aj'(0) and obtain 



,2 An) 



,2c(n) 



(34b) 



(35) 



In order to determine the other two eigenvalues A ; 1 
and A^ 2 , which depart from zero with Og, we set flp = 



^Note that unlike t he off-resonance case with S a 2 = Sp2 = 
S a p = (Sec. Ill A 4), no additional zero eigenvalues exist for 



A„ = when S$ = 5$ = 0. 







in Eq. (|22|) and divide by A (thus removing the root A~ ) . 
Keeping the terms of lowest order with respect to fig and 
A, we find that — r2|/3^ + A 2 « 0, and we identify A^ and 
Ar 2 as the two roots of this equation, 



(36) 



2. Late-time eigenvalues 

In a similar fashion, we find that for t — > +oo, the 
three vanishing eigenvalues behave as 



A+ 



1 



a 2 n Sp - 2a n p n S { $ + fiS$>] Q% (37) 



Vi ~ " a n^p> A+2 « a n tt P . (38) 



3. Connectivity 

It is straightforward to show that the eigenstate associ- 
ated with A J tends to state ipi as t — ► — oo and the eigen- 
state associated with A+ tends to state ipf as t — > +oo. 
The eigenstates corresponding to the "large" eigenvalues 
tend to superpositions of states ipf and ip n initially and 
to superpositions of states ipi and ip n finally. Hence, if 
A~ and A+ correspond to the same eigenvalue, the corre- 
sponding eigenstate will be the desired AT state <£>t(*)- 
We have seen above that in the general off-resonance case, 
this may or may not take place. In the present case of a 
single-photon resonance, however, this is always the case. 
To show this we first note that the eigenvalues, which do 
not vanish at ±oo, do not interfere in the linkages be- 
tween the vanishing eigenvalues because each of the non- 
vanishing eigenvalues Xk(t) tends to the corresponding 
detuning at both ±oo. Hence, the eigenvalues that 
are above (below) the three vanishing eigenvalues at — oo 
remain above (below) them at +oo as well. 

Let us now consider the linkages between the three 
vanishing eigenvalues. Insofar as £lp/Qs —* as t — > 



This means 



oo, we have X t , < A s < X t 2 . Also, since Clg/Clp 



as t 



we have A, , < At < A 



1.2- 



A u , A; 



Xf , and A 



1.2 



, 2 take place, and therefore, the AT state frit) always 



that the linkages A ; 1 
A' 

exists when the lasers are tuned to resonance with an 
intermediate state, which is indeed seen in Figs. |^ and]?] 



C. Examples 

In Fig. U we have plotted the time evolutions of the 
three eigenvalues that vanish at ±oo (upper row of fig- 
ures) in the case of N = 2 intermediate states for 
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FIG. 8. The upper row of figures show the time evolu- 
tions of the three eigenvalues that vanish at ±oo in the case 
of N — 2 intermediate states for two combinations of coupling 
strengths: cti = /3i = 1, a?2 = 02 — 0.5 (upper left figure) and 
qi = Pi — 1, Qf2 = 2, and (3% = 0.5 (upper right figure). 
The detunings are the same in both cases: Ai = 0, A2 = Oo. 
The solid curves are calculated numerically and the dashed 
curves show our asymptotic approximations (pB[)-(|3^). The 
bottom row of figures show the components of this eigenstate 
which is equal to the bare state ipi initially. It corresponds 
to the eigenvalue which is shown by a thick curve in the cor- 
responding upper figure. This is the zero eigenvalue Xd for 
the left figure and the eigenvalue At whose asymptotics at 
early times is described by AJ for the right figure. The labels 
denote the bare states to which the components belong and 
the final-state component is shown by a thick curve. 



two sets of coupling strengths and the same detunings, 
Ai = 0, A2 = Sloj he., the lower intermediate state is 
on single-photon resonance. The top left figure is for 
proportional coupling strengths (|l8| ) which give rise to a 
zero eigenvalue Xd and correspondingly, to a dark state 
tfiD(t). The top right figure is for a case when Eq. jl^ ) 
is not satisfied and there is no zero eigenvalue. Note the 
pair of eigenvalues A^ and \J 2 which depart from zero 
in opposite directions with the arrival of the Stokes pulse 
at early times and the pair of eigenvalues A^ and A ; + 2 
which vanish with the disappearence of the pump pulse 
at late times. The bottom row of figures show the com- 
ponents of this eigenstate which is equal to the bare state 
ipi initially; it corresponds to the zero eigenvalue for the 
left figure and to the eigenvalue whose asymptotics at 
early times is described by Xj for the right figure. The 
squared components of this eigenstate give the popula- 
tions of the four bare states in the adiabatic limit. As 
predicted by our analysis, the AT state is seen to exist in 
both cases. The difference is that for the left column of 
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FIG. 9. The final-state population Pf as a function of 
the pulse width T in the case of TV = 2 intermediate states 
for the same two sets of interaction parameters as in Fig. [| 
The solid and dashed curves correspond to the left and right 
columns in Fig. ^, respectively. The detunings are the same in 
both cases: Ai = 0, A2 = S7o- The solid curve is for coupling 
strengths on = /?i = 1, 02 = #2 = 0.5 (then the AT eigenvalue 
is Ad = 0), and the dashed curve is for ai = /3% = 1, Q2 = 2, 
and P2 = 0.5 (then the AT eigenvalue is nonzero). 



figures, the AT state is the dark state (||), which does not 
contain components from the intermediate states, while 
for the right column of figures, the AT state contains such 
components. 

In Fig. ||, we have plotted the final-state population Pf 
as a function of the pulse width T in the case of N = 2 
intermediate states for the same two sets of interaction 
parameters as in Fig. g. The solid and dashed curves in 
Fig. H correspond to the left and right columns in Fig. [|, 
respectively. In agreement with our conclusions, an AT 
state exists in both cases and the final-state population 
Pf approaches unity as T increases. 



V. DEGENERATE RESONANT INTERMEDIATE 
STATES 

Let us suppose now that No detunings vanish, i.e., that 
there are No degenerate resonant intermediate states, 
and let us assume without loss of generality that these 
states are tp± , ip2 > ■ ■ ■ , fpN ■ ^ two °f the detunings are 
equal to zero, Ai = A2 = 0, we have 



dctH = fi£n!x>i 2 (ai/?2 - a 2 Pi) 2 , 



(39) 



and hence, a zero eigenvalue exists when a\/ (3\ = ail f$i- 
If three or more detunings are equal to zero then det H = 
0, and a zero eigenvalue always exists, with no restrictions 
on the interaction parameters. 



A. Proportional couplings 

In the zero-eigenvalue eigenstates, the components 
di(t) and df{t) of states ipi and ipf satisfy the equations 



a n n P {t)a t {t) + /3 n n s (t)a f (t) = 0, 



(40) 



with n — 1,2, ... , No- A nonzero solution for di(t) and 
a f (t) requires that 



«jVo 
Pn 



= 1. 



(41) 



Otherwise, a zero-eigenvalue eigenstate cannot be an AT 
state. If relation ( |4l| ) is fulfilled, the amplitudes of the 
degenerate states are linearly dependent]^] 



Cn(t) = Ci(t)- 



(n = 2,3,...,/Vo). 



(42) 



This relation allows to replace in the Schrdinger equation 
(|l|) the probability amplitudes of the degenerate states by 
an effective amplitude given by 

Ceff(t) =/iCi(t), (43) 

with pump and Stokes Rabi frequencies given by 

fij»eff(t) = MM*)> ^S.effO) = ^S(t), (44) 



where /1 = [\/ai)Ja\ 



. Thus the orig- 



inal problem with Nq resonant states is reduced to an 
equivalent problem involving a single resonant state. As 
we pointed out in Sec. an AT state <pr(t) always ex- 
ists in this case. Moreover, in the case when all couplings 
(and not only those for the degenerate states) are pro- 
portional the AT state is the dark state ipD(t). 



B. Arbitrary couplings 

If Eq. @ is not fulfilled, the zero-eigenvalue eigen- 
state^) cannot be an AT state because the components 
a,i(i) and a,f(t) from states ipi and ipf vanish. There still 
might be a possibility that one of the two nonzero eigen- 
values, which vanish at ±00, corresponds to an AT state. 
We shall show, however, that this is not the case. 

For No degenerate resonant states, it can readily be 
shown that the number of zero eigenvalues is No — 2 for 
Q P (t) ^ and Q s {t) ^ 0, N for Q P (t) = and Q s (t) ^ 
[or for fip(i) ^ and tt s (t) = 0], and A^ + 2 for 
Qp(t) = £ls(t) = 0. This means that when the Stokes 
pulse arrives at early times, it lifts the degeneracy of 



§ This is so because the function f(t) = a n C\(t) — aic n (t) 
satisfies the differential equation df(t)/dt — with the initial 
condition /( — 00) = 0; hence, f(t) — 0. 
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two of the A*o + 2 zero eigenvalues. When the pump 
pulse arrives later, if lifts the degeneracy of another pair 
of the remaining No eigenvalues. The reverse process 
occurs at large positive times. The remaining Nq — 2 zero 
eigenvalues stay degenerate all the time. The implication 
is that the initial state ipi and the final state ip / are given 
by superpositions of eigenstates both at — oo and +00, 
which means that an AT does not exist. 



VI. ADIABATIC ELIMINATION OF THE 
OFF-RESONANCE STATES 

A. The off-resonance case 

An insight into the population transfer process can be 
obtained from the adiabatic-elimination approximation. 
When the single-photon detuning of a given interme- 
diate state ipk is large compared to the couplings f2p,fe 
and fis.fc of this state to states ipi and ipf, this state 
can be eliminated adiabatically by setting dck/dt = 
and expressing Ck from the resulting algebraic equation. 
By adiabatically eliminating all intermediate states the 
general (N + 2)-state problem is reduced to an effective 
two-state problem for the initial and final states, 
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FIG. 10. The final-state population Pf as a function of 
the pulse width T in the case of N = 3 intermediate states for 
three combinations of detunings, denoted in the figure. The 
upper T-scale is for the dotted curve while the lower T-scale 
is for the dashed and solid curves. The coupling strengths 
are the same in all cases, «i = /3i = 1, Q2 = 0.6, 02 = 1.0, 
a 3 = 1.2, #, = 0.6. 
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(45) where t c is the crossing point. For the Gaussian shapes 
d) we have t c = (T 2 /8r) ln(5« 2 /S a 2 ) and 



incorrect and that an AT exists even in this case, as long 
as S a 2S/32 > 0. This somewhat surprising discrepancy 
derives from the fact that for S a p = 0, the effective cou- 
pling between ipi and ipf is so small that it is lost in the 
course of the approximation. Hence, this approximation 
provides a useful hint for the least favorable combina- 
tion of parameters which results in the weakest effective 
coupling between ipi and ipf, consequently, the adiabatic 
limit is approached most slowly in this case. 

The adiabatic-elimination approximation allows to es- 
timate how quickly the adiabatic limit is approached 
when the AT state exists. Then, as we noted above, we 
have a level-crossing transition, the probability for which 
can be roughly described by the Landau-Zener formula, 



with 



The "detuning" in this two-state problem is A e g(i) = 
(t)Sp2 — Q p (t)S a 2. Obviously, if S^S^ > 0, A e g(t) 
has different signs at ±00 and the transition is of level- 
crossing type, while if S^S^ < 0, A B g(t) has the same 
sign at ±00 and there is no crossing. Hence, in the 
adiabatic limit, the transition probability from state ipi 
to state ipf will be unity for S a 2Sp2 > and zero for £ 
S a 2Sp2 < 0, in agreement with the AT condition (pl|). 

The "coupling" in the effective two-state problem ( Ma ) 
is Q e ff(t) — Hp(t)fls(t)S a 0. Obviously, it vanishes for 
Sa/3 = which suggests that there is no transition from 
state ipi to state -0/, both for S a 2Sg2 > and S a 2Sp2 < 0. 
However, we know from Sec. HI A 4 that this prediction is 



"c 2 ff(*c) 

A cS (t c ) 



(tt T) 2 £, 



T 



<? 2 



At 



\J S a 2 S t 



: exp 



2 



2T 2 
J^2 



T 2 
32t~ 2 



In- 



(47) 



(48) 



The larger this parameter, the faster the adiabatic limit is 
approached. We thus conclude that from the adiabaticity 
viewpoint, the most favorable case is when S a 2 = S(j2 
and the ratio S^p/^S^S^ is large. Not surprisingly, 
the Landau-Zener parameter ([47j ) is also proportional to 
(f^o? 1 ) 2 , which is essentially the squared pulse area. 

In Fig. [To] the final-state population Pf is plotted 
against the pulse width T in the case of N = 3 inter- 
mediate states for three combinations of detunings. The 
coupling strengths are the same in all cases. The parame- 
ters for the dotted curve are chosen so that S a p = while 
those for the other curves ensure that S a p 7^ 0. The pa- 
rameter ( fl8| ) is £ = for the dotted curve, £ ~ 0.326 for 
the dashed curve, and £ « 1.367 for the solid curve. As a 
result, the adiabatic limit is approached most slowly for 
the dotted curve, more quickly for the dashed curve, and 
most quickly for the solid curve. 



P 



f 



1 



_ e -^cff(tc)/Aeff(tc) 



(46) 
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B. The on-resonance case 

When a certain intermediate state ip n is on single- 
photon resonance, A„ = 0, it cannot be eliminated adi- 
abatically. By eliminating all other intermediate states, 
the general TV +2- state problem is reduced to an effective 
three-state problem, 



dt 



' Ci ' 








- c f- 





q2 c(") 



a n £lp 













- c f- 



(49) 



where the sums are defined by Eqs. (|3l|). Compari- 
son with the standard three-state STIRAP shows that 
the off-resonant states induce dynamic Stark shifts of 
states ipi and ipf, which result in a nonzero two-photon 
detuning. Moreover, the off-resonant states induce a di- 
rect coupling between states ipi and ipf. Careful exami- 
nation of Eq. ( pgj ) shows that the AT state ( p~5| ) always ex- 
ists, but it involves in general a nonzero contribution from 
the intermediate state ip n . This contribution vanishes 
when the proportionality condition ( |l8| ) is fulfilled; then 

— iSfc' = S„) and the AT state is the dark state, 

a z p z otfj ' 



<PT(t) = tfD{t)- In particular, if = S 
the multi-A system behaves exactly like STIRAP. 



,(n) 

2 



;(n) 
3 a 



0, 



VII. DISCUSSION AND CONCLUSIONS 

We have presented an analytic study, supported by nu- 
merical examples, of adiabatic population transfer from 
an initial state ipi to a final state ipf via TV intermediate 
states by means of two delayed and counterintuitively or- 
dered laser pulses. Thus this paper generalizes the origi- 
nal STIRAP, operating in a single three-state A-system, 
to a multistate system involving N parallel A-transitions. 
The analysis has shown that the dark state if nit), which 
is a linear combination of ipi and ipt and transfers the 
population between them in STIRAP, remains a zero- 
eigenvalue eigenstate of the Hamiltonian (^) only when 
condition ( jlq ) is fulfilled. Hence, in this case the multi-A 
system behaves very similarly to the single A-system in 
STIRAP. Condition (fl8|), which is essentially a relation 
between the transition dipole moments, requires that for 
each intermediate state ipk, the ratio £lp,k(t)/£ls,k{t) be- 
tween the couplings to states ipi and ipf is the same and 
does not depend on k. Moreover, this condition ensures 
the existence of the dark state both in the case when all 
intermediate states are off single-photon resonance and 
when one or more states are on resonance. 

When condition ( |l8| ) is not fulfilled the dark state 
ipoit) does not exist but a more general adiabatic- 
transfer state fxit), which links adiabatically the ini- 
tial and final states ipi and ipf, may exist under certain 



conditions. Unlike ipz>(t), state frit) contains contri- 
butions from the intermediate states which therefore ac- 
quire transient populations during the transfer. We have 
shown that when one and only one intermediate state is 
on resonance, the AT state always exists. When more 
than one intermediate states are on resonance, the AT 
state exists only when the proportionality relation flis| ) 
is fulfilled, at least for the degenerate states. In the off- 
resonance case, the condition for existence of ipr(t) is 
given by Eq. ( p9|) which is a condition on the single- 
photon detunings and the relative coupling strengths. 
It follows from this condition that when the pump and 
Stokes frequencies are scanned across the intermediate 
states (while maintaining the two-photon resonance), the 
final-state population Pf passes through N regions of 
high transfer efficiency (unity in the adiabatic limit) and 
N — 1 regions of low efficiency (zero in the adiabatic 
limit). Each of the low-efficiency regions is situated be- 
tween two adjacent intermediate states, while each in- 
termediate state is within a region of high efficiency, as 
shown in Figs. |] and 0. Our results suggest that it is most 
appropriate to tune the pump and Stokes lasers either 
just below or just above all intermediate states because 
there firstly, the AT state always exists; secondly, the 
adiabatic regime is achieved more quickly; thirdly, the 
transfer is more robust against variations in the laser pa- 
rameters; and fourthly, the transient intcrmedaite-state 
populations, which are proportional to A^ 2 , can easily 
be suppressed. 
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